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We present an update on our previous determination of the Lambda parameter in QCD. The main emphasis 
is on results for two flavours of light dynamical quarks, where we can now almost double the amount of data 
used, including values at smaller lattice spacings. The calculations are performed using O(a) improved Wilson 
fermions. Little change is found to previous numerical values. 



The A parameter is one of the fundamantal pa- 
rameters of QCD, setting the scale for the run- 
ning coupling constant a s . In this contribution 
we shall update our previous work, P^, both for 
quenched (nt = 0) and unquenched (nt = 2) 
0(a) improved Wilson ('clover') fermions. Specif- 
ically we are now able to use for 

• quenched fermions, the force scale ro/a up 
to P = 6.92, (previously (3 < 6.4), 

• unquenched fermions, improved statistics 
and additional quark masses at the previous 
(3 values of 5.20, 5.25, 5.29 for ro/a together 
with new results at j3 = 5.40 (at three quark 
masses). 

The 'running' of the QCD coupling constant as 
the scale changes is controlled by the /3-function 



renormalisation having introduced a scale M to- 
gether with a scheme S. Integrating this equation 
gives 
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where, perturbatively 

P s (gs) = -b g 3 s - b l9 % 



b S 2 gl - bfg% 



'Talk given by R. Horsley at Lat04, Fermilab, USA. 



where A 5 , the integration constant, is the funda- 
mental scheme dependent QCD parameter. Re- 
sults are usually given in the MS scheme, with 
the scale M being denoted by In this scheme 
the first four /3-function coefficients are known, 
bf s being found in The running coupling 

a^' s (fi) = grrsifJ-) 2 1 '4tt is plotted in Fig. ^ for 
nf = 2, by solving the previous equation (nu- 
merically) using successively more and more co- 
efficients of the /3-function. The figure shows 
an apparently fast convergent series (cf 3- to 4- 
loop), certainly in the range we are interested in, 
n/A MS ~ 8. A very similar result holds for ?!/ = 
but with slightly lower curves. 
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Figure 1. af s (fi) versus fi/A MS for nf = 2. 

On the lattice we also have a A parameter, 
a A D -F D (.g D (a)), 

where to help convergence of lattice perturbative 
expansions we use = g 2 (a)/uQ with Uq the av- 
erage plaquette value. To calculate A MS , we shall 
compute gjfs at some appropriate scale /i* from 
gn (a) and then using the tq scale, extrapolate 

r A~ ee (^) F^fe(//))a// , 

to the continuum limit. 

Equating lattice and continuum expressions 

[F^igjrs^W 1 = a^lF^g^a))}- 1 , 

and expanding as 
1 1 

9a(a) 

[2bo In afi — t°] + [2b i In a[i — i° ]<?g (a) . . . , 

gives = 26 In A~/A D and b° = 6p + M? - 
60^2 • For (hopefully) good convergence of this 
series we choose the scale so that the O(l) term 
vanishes, afi* = exp(t°/26o)- 

For ti the general expression is known for nf, 
c sw and linear terms in n/am g , while for the 
nfam q dependence is not known, and refer- 
ences therein. We can estimate the scales as 
ft* = 2.63/a, nf = and ft* ~ 1.4/a for u/ = 2. 
£3 (the ctes(/i) 4 , hi afi independent term) is not 



known. So equivalently 63 is not known. However 
a Pade estimate gives bf m (frf) 2 /^, and is small 
and in reasonable agreement with the known co- 
efficient in the MS scheme, £Q. Assuming this 
also holds for 63 gives little change to the results 
presented here. For complete 0(a) cancellation, 
[2], we need g 2 — g 2 (l + b g am q ) where pertur- 
batively b g = 0.01200n /ff 2 + 0{g A ), which with 
c S w — 1 + 0(g 2 ) then gives no mass dependence 
in t±. This indicates little quark mass depen- 
dence in the fit formulae (indeed there is more in 
the numerical data). Finally to further improve 
the convergence of the series, we tadpole improve 
the tf coefficients — c sw Uq (for t° + t^g'u) 
further reducing the size of the n / term in t% ■ 

In Fig. |2 we show the quenched {rif =0) re- 
sults. The data lies on a straight line (as a func- 
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Figure 2. roA MS versus (a/ro) 2 for rif = 0, together 
with a linear extrapolation to a = 0. The last point has 
not been included in the fit. 

tion of (a/ro) 2 ) at least over a -1 ~ 2 - 6.5 GeV 
or ft ~ 5 - 17 GeV, using the value for r of 
ro = 0.5 fm. This gives a result of roA MS = 
0.613(2)(25) or A"^(0) = 242(1)(10) MeV where 
the first error is statistical and to estimate the 
systematic uncertainty, the second error takes a 
g 4 coeff. = 25% x g 2 coeff. (which is very much 
greater than when using the Pade 63 estimate). 

For unquenched (rif — 2) fermions, due to 
the sea quark, the fit ansatz is not so sim- 
ple as we must consider both chiral and con- 
tinuum extrapolations. We take for finite a, 
aA MS | mg ^ ,a#o = aA MS \ mq=0 ^ + Dam q + . . . or 
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n)A MS | m ^ ,a#o = r A MS | m<!= o,^o+from,+. . .. 
After chiral extrapolation we would thus expect 

r A~\„ lq = ^o = r A™\ m ,=Q,a=Q + B(a/ P ) 2 + 
. . . with p = r \ m - . Together with (a/r ) 2 = 



(a/p) + Eam q 



this gives our fit ansatz as 



r A MS = A + B(a/r ) 2 + Cam q + Dr m q . 

So by subtracting out the B and C terms from 
roA MS we can consider the chiral extrapolation 
and similarly by subtracting out the C and D 
terms we may consider the continuum extrapola- 
tion 1 . In Fig. |3 we show the results, a -1 ranges 
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Figure 3. roA MS versus rom q (upper picture) and ver- 
sus (a/ro) 2 (lower picture) for rif = 2, together with ap- 
propriate extrapolations (am q from Wt). 

at least over a -1 ^2-3 GeV or fi ~ 3 - 4 GeV 
This gives a result of r A~ = 0.563(10) (70) or 
A~(2) = 222(4)(28)MeV where again the first 
error is statistical and the second error is obtained 
by taking a g A cocff. = 25% x g 2 coeff. which again 
is much larger than the error found when using a 
Pade 63 estimate, setting c sw = 1 + 0(g 2 ) or in- 
cluding an additional (am q ) 2 fit term. Note that 
this result is consistent with that obtained in .6; . 

Finally in Fig. 0] we present results for different 
n/ . Our result lies somewhat low in comparison 
with phcnomcnological results. Alternatively us- 
ing the matching procedure as in we find for 
71/ = 5, af(m z ) = 0.1084(6)(38).' ' 

1 An alternative procedure is first to extrapolate both rrj/a 
and Uq to the chiral limit, evalute roA MS and then extrap- 
olate to the continuum limit; this gives similar results, 
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Figure 4. A MS (n^) versus nf The open diamonds are 
from using af^{M z ) = 0.1183(27) (n f = 5) to match 
to nf = 4 and nf = 3, while the open square is from 
The filled circles are the results reported here. 
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